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The measurements of very low level signals at 
low frequency is a very difficult problem, because 
environmental noise increases in this frequency do- 
main and it is very difficult to filter it efficiently. 
In order to counteract these major problems, we 
propose a simple and generic denoising technique, 
which mixes several features of traditional feedback 
techniques and those of noise estimators. As an ex- 
ample of application, large band measurements of 
the thermal fluctuations of a mechanical oscillator 
are presented. These measurements show that the 
proposed denoising technique is easy to implement 
and gives good results. 

1 Introduction 

The measurements of very low level signals at very 
low frequency (VLF) is a very difficult problem, be- 
cause environmental and electric / magnetic noises 
often increase in this frequency domain. Further- 
more, it is well known that the difficulties of iso- 
lating an experimental setup from these unwanted 
noise sources increase by reducing the measuring 
frequency. Typical examples are the screening of 
low frequency magnetic fields or the isolation of a 
measurement from the unwanted environmental vi- 
brations. Many techniques have been proposed and 
accurately applied to reduce the effects of the un- 
wanted noise sources. 

The simplest techniques are of course the passive 
ones. Let us consider in some details the problem 
of vibration isolation (the magnetic field screening 
presents similar problems). 

In a typical laboratory environment, vibrations 
transmitted through the floor normally have a fre- 
quency spectrum from dc to a few hundred Hz. The 
reduction of these vibrations is usually obtained by 



installing the experiment on floating tables which 
have horizontal and vertical resonance frequencies 
around f Hz. Thus, noise reduction is obtained 
only for frequencies larger than the natural reso- 
nant frequencies of the table. This is, of course, 
an excellent method for high frequency measure- 
ments, but at frequencies close to and smaller than 
f Hz this method becomes useless. Exactly at res- 
onance, noise is even enhanced. To overcome these 
problems that appear at VLF, feedback techniques 
have been used. These techniques require detec- 
tors which measure the noise signals and actuators 
which reduce the acceleration of the table plate. 
Similar techniques are of course used for screen- 
ing VLF magnetic fields. Indeed, high frequency 
magnetic fields are screened by Faraday cages and 
the VLF components are subtracted by a feedback 
technique [l]. 

Feedback techniques are widely used, but they 
are limited by the noise of the detectors and of 
the actuators. Their calibration is often very com- 
plex and requires very tedious operations in order 
to work properly. Furthermore, they can be only 
applied in all of the cases where the environmen- 
tal noise can be accurately measured. When this is 
not possible and / or one is not interested in stabi- 
lizing a system on a given working point, but only 
in reducing the noise on a given signal, other tech- 
niques may be used. From the signal analysis point 
of view, most of the techniques that have been pro- 
posed and accurately applied to reduce the effects 
of the unwanted noise sources rely upon the knowl- 
edge of the response function of the system under 
study, and a guess about the noise, which is of- 
ten supposed to be a random variable belonging to 
some class of signals, e.g. ergodic and second order 
stationary signals El El El El- However, this guess 



is often limited. 

The simple denoising technique proposed in this 
paper actually combines several aspects of the two 
above mentioned techniques. Indeed, in many cases 
one has access to the environmental noise, but one 
does not need to stabilize a setup on a given work- 
ing point, but only to reduce the effect of the noise 
on a given signal. Therefore, rather than carrying 
out a somewhat sophisticated and expensive feed- 
back system, which could even fail to solve the noise 
problem at VLF, we developped a simple denoising 
technique whose principle precisely lies in measur- 
ing the residual noise when passive "screening" de- 
vices are already used. The general principles and 
the limits of the new technique are described in 

SecEl 

This denoising technique was motivated by the 
study of the violation of the fluctuation dissipation 
theorem (FDT) in out of equilibrium systems as 
aging glasses. This is a subject of current interest 
which begins to be widely studied in many differ- 
ent systems G3 El EH EI] Therefore, in Sec.E|we 
propose an application of this new technique to the 
measurement of VLF mechanical thermal fluctua- 
tions. The experimental results, presented in Sec.QI 
show the quality of the noise reduction. Finally, in 
Sec.0 we discuss some other possible applications 
and we conclude. 

2 A simple denoising tech- 
nique 

Suppose one has to measure a signal on which an 
external noise is superimposed. Let us call this sig- 
nal the true signal Xtrue(t) and the external noise 
x C nv(t), so that the measured signal can be written 
in the additive manner 



x(t) = x tIUC (t) + x cm (t). 



(1) 



The only assumptions needed are that x tIU e(t) and 
%env(t) are uncorrelated and both stationary pro- 
cesses (as we will see, the hypothesis of stationarity 
can be weakened), and can be written in an addi- 
tive manner like in Eq.^ In addition, we assume 
that one can directly measure x, whereas x cnv is 
measured with a detector whose output signal is 
Xdet(t). If the noise is small, in the limit of lin- 
ear response theory, it can be stated that Xdct is 



linearly related to x cnv , which means that there ex- 
ists a hypothetical response function -R~j~ (w) such 
that 

£det(w) = Rd 1 (w)£onv(w), (2) 

where f(uj) = J R f(t)e~ ZUJt dt is the Fourier trans- 
form of f(t). Notice that no hypothesis is done on 
Rd, which is in principle unknown. 

The Fourier transform of the true signal thus 
reads 

Xtruc = X- i?d£dct- (3) 

Assuming that the true and external noise signals 
are uncorrelated, that is (xtruc ^dct) — Oi one can 
compute the kernel Rd as 



Rd 



[xx. 



dct/ 



Fdct 



(4) 



where (•) stands for the ensemble average. Thus 
Eqs.El and fallow us to compute the signal and its 
spectrum: 



Ztrue(*) = X(t) 



fldMid-tHe*"*^, (5) 



and 



(\xtn 



(1*1 



|-Rd| 2 (|^dct| 2 ) 



(6) 



Therefore, Xtrue can be computed from the simul- 
taneous measurements of x and x cnv . 

We see that the hypothesis of stationarity is not 
really necessary, because x, R^ and Xdet can be 
slowly varying functions of time, with a character- 
istic time t. In such a case, if the ensemble average 
is performed in a time T such that T <C t, then 
Eqs.ElEEland|Hlcan be still applied on intervals of 
length T. This observation makes this simple tech- 
nique very powerful, because the response of the 
system to the environmental noise can change as a 
function of time and of the external noise source. 
Thus, Rd is a dynamical variable which can be com- 
puted in each time interval of length T, and which 
allows to retrieve the true signal. 

However, the signal x rcc reconstructed using 
Eqs.03 EJ and HI will differ from Xtruc because of 
experimental errors. One source of error is the 
noise of the detectors and of the amplifiers, which 
introduces an extra additive noise term r)(t) in 
Eq.[IJ which is uncorrelated with Xtrue and Xdot, 
thus x YCC (t) = x tIUC (t) + i](t). However, r\ can be 



done very small and it does not constitute the main 
source of error. The main one is the limited num- 
ber N of ensemble averages that can be done in the 
time T. This is very important because if Rd is a 
slowly varying function oft, then one has to impose 
T -c t in order to retrieve the true signal. Finally, 
it has to be pointed out that the advantage of the 
technique is when the amplitudes of Xt ruo and Xdet 
are comparable, that is, when the signal to environ- 
mental noise ratio is either smaller than or equal to 
one. 





The accuracy of the technique has been checked 
on several artificial signals. We have chosen for 
xt r uc and Xdot either random or periodic signals. 
The random signals may be either colored or white 
noise with Gaussian or uniform distribution. To 
estimate the error of the reconstruction, we first 
consider the difference A between the reconstructed 
x Iec and xtrue, that is A(t) = x ICC (t) -x tTUC (t). We 
then compute the ratio i? rec between the rms of A 
and that of Xt ruc , which is a good indicator of the 
error of the reconstructed signal. As expected we 
find that R lcc ~ 1/y/N for large N. 



Figure 1: Artificial signal, (a) True signal as 
a function of time. In this example Xtme(t) — 
sin |ji + sin Jgjjt. (b) The signal iasa function of 
time. The noise is a Gaussian white noise of vari- 
ance 4. (c) The reconstructed signal as a function 
of time, (d) Residual error A of the reconstruction 
as a function of time. 



Application to thermal fluc- 
tuations measurements 



An example of the reconstruction is given in 
Fig.Q] We see that although the signal is com- 
pletely erased by the noise (cf Fig. lb) the recon- 
struction is quite good. It is obvious that this is an 
extremely simple example, but as we will see the 
technique becomes very interesting when Rd is a 
slowly varying function of t. 

To conclude this section, it should be stressed 
that, from the signal analysis point of view, we 
derived a method in a way similar to the Wiener 
filtering, which aim at separating (in an optimal 
sense, see El E3) two random signals x trU c an 
Xdet, which are supposed to be ergodic second order 
stationary and uncorrelated random signals, and 
can be written in an additive manner like in Eq.Q] 
Then, we extended this denoising technique to non- 
stationary signals in a simple and original man- 
ner. In a more general study, nonstationary signals 
could be addressed to the Kalman filtering (also re- 
ferred as to the Kalman-Bucy filtering) , which can 
be considered as the extension of the Wiener filter- 
ing to the case of nonstationary signals ^1 El IE]- 



In this section, we describe a useful application of 
this technique to the measurement of thermal fluc- 
tuations of a mechanical oscillator, whose damp- 
ing is given by the viscoelasticity of an aging poly- 
mer glass. This is an important experimental mea- 
surement which is extremely useful in the study of 
the violation of the FDT in out of equilibrium sys- 
tems, specifically in aging glasses. This violation is 
a subject of current interest which has been stud- 
ied mainly theoretically |1-'U I14j . However, there 
are not clear experimental tests of these theoretical 
predictions, which have to be checked on real sys- 
tems by studying the VLF spectra of mechanical 
thermal fluctuations. Thus, the main purpose of 
our study is to have a reliable measurement of this 
VLF spectra in an aging polymer. 

To study this spectrum, we have chosen to mea- 
sure the thermally excited vibrations of a plate 
made of an aging polymer such as Polycarbonate. 
The physical object of our interest is a small plate 
with one end clamped and the other free, i.e. a can- 
tilever. The plate is of length I, width a, thickness 
b, mass mp i yc . On the free end of the cantilever 
a small golden mirror of mass m m i rr0 r is glued. As 



described in the next section, this mirror is used 
to detect the amplitude x c of the transverse vibra- 
tions of the cantilever free end. The motion of the 
cantilever free end can be assimilated to that of a 
driven harmonic oscillator, which is damped only 
by the viscoelasticity of the polymer. Therefore, 
the Fourier-transformed equation of motion of the 
cantilever free end reads 



-mui + K(lu)]x c = F c: 



(7) 



where x c is the Fourier transform of x c , rn is the 
total effective mass of the plate plus the mirror, 
K = K' + iK" is the complex elastic stiffness of the 
plate free end, and F ex t is the Fourier transform 
of the external driving force. The complex K(uS) 
takes into account the viscoelastic nature of the 
cantilever. From the theory of elasticity Q2] ° ne 
obtains that, for VLF, excellent approximations for 
to and K are: 



and K 



(3.52) 2 
Eab 3 



^Polyc i "^mirror; 



4/ 3 



(8) 
(9) 



where E = E' + iE" is the plate Young modulus. 
Notice that if m m i rror = 0, then one recovers the 
smallest resonant frequency of the cantilever [T5] . 
For Polycarbonate at room temperature, E is such 
that E' = 2.2 x 10 9 Pa and E" = 2 x 10 7 Pa, and 
its frequency dependence may be neglected in the 
range of frequency of our interest, that is from 0.1 
to 100 Hz [T^- Thus we neglect the frequency de- 
pendence of K in this specific example. 

When F cxt = 0, the amplitude of the thermal 
vibrations of the plate free end xt is linked to its 
response function x vl & the FDT [Tfi| : 



2 2k B T A 
(Ft| ) = Imx, 



(10) 



where (\xt\ ) is the thermal fluctuation spectral 
density of x c , k B the Boltzmann constant and T 
the temperature. From Eq.0one obtains that the 
response function of the harmonic oscillator is 



1 



-Fcxt m[uj 2 -uj 2 -i (sign w) juj 2 
where lu 2 = K'/m and 7 = K"/K'. 



(11) 



Inserting Eq.^Jinto En . 1101 one can compute the 
thermal fluctuation spectral density of the Polycar- 
bonate cantilever for positive frequencies: 



(Ft I ) = 



7^o 



uj m[(ujQ 2 — uj 2 ) 2 + (7W0 2 ) 2 ] 



(12) 



Notice that (\£t\ ) ~ ^ _1 for uj <C Wo, because 
the viscoelastic damping K" is constant in our fre- 
quency range. In the case of a viscous damping (for 
example, a cantilever immersed in a viscous fluid) 
K" = aoj, where a is proportional to the fluid vis- 
cosity and to a geometry dependent factor. Then 
the spectrum of the thermal fluctuations of the can- 
tilever free end, in the case of viscous damping, is 



(\xt\ 



2k B Ta 



m 2[( W g-W a ) 2 + (» S 



(13) 



which is constant for uj <C ojq. Therefore the fluc- 
tuation spectrum shape depends on K"(uj). In the 
case of a viscoelastic damping (see Eq.El, the ther- 
mal noise increases for uj <C uj , and with a suitable 
choice of the parameters the VLF spectrum of an 
aging polymer can be computed using this method. 
However, the cantilever is also sensitive to the 
mechanical noise, and the total displacement x c of 
the cantilever free end actually reads x c = xt + 
x acc , where x acc is the displacement induced by the 
external mechanical noise. Thus, it is important 
to compute the signal-to-noise ratio of our physical 
apparatus, which we define as the ratio between 
the thermal fluctuations and the mechanical noise 
spectra. To compute the latter, we consider that 
the support of the cantilever is submitted to an 
external acceleration a cxt , whose Fourier transform 
is a ex t. We rewrite Eq.^lwith F e xt = rna-ext, which 
yields 

Ocxt 



UJ Z 



UJ^ 



1 ' 



(14) 



where x acc is the Fourier transform of x acc . Far 
from the resonance frequency, that is for uj <§; 
ujo, one has x acc ~ a C xt/wo 2 , which finally yields 
l^accl ~ |flcxt| /^o 4 j whereas the thermal fluc- 
tuation spectral density of x reads (\&t\ ) ~ 
=--^- Therefore, the signal-to-noise ratio 



reads 



(\x T h 2k B T juj 2 



v TO(|a oxt | 2 )' 



(15) 



which is proportional to 



-yEab* 



m 2 i3(|« acc j 2 > 



, for u> <C U)q. 



Towards the Nomarski Interferometer 
Fixed Mirror 



Notice that the signal-to-noise ratio of Eq.^1 in- 
creases if the set of parameters {a, b, I, m} is opti- 
mized to make uj as large as possible within the 
frequency range of interest, and within the experi- 
mental constraints. 

Let us estimate the amplitude of y (\xt\ ) at 
v = u>/2ir = 1 Hz for the following choice of the 
parameters: 7 ~ 10 -2 , I ~ 10 mm, a ~ 1 mm, 
b = 125 /im and m m i rror < 1CT 3 g. We find 

v Q ~ 100 Hz and y 1 '{\x T (l Hz)| 2 ) ~ 10~ n m/VSz, 
which is a very small signal. As a consequence, ex- 
tremely small vibrations of the environment may 
greatly perturb the measurement. Therefore, to 
increase the signal-to-noise ratio of the measure- 
ments, one has to reduce the coupling of the can- 
tilever to the environmental noise (acoustic and 
seismic) using vibration isolation systems. This 
may be not enough in this specific case because 
of the the smallness of the thermal fluctuations. 
Then we have applied the technique described in 
the previous section in order to recover xt from 
the measurement of x c . The experimental results 
are described in the next section. 



4 Experimental results 

The measurement of x c is done using a Nomarski 
interferometer (for detailed reviews, see |18M19ll20] l 
which uses the mirror glued on the Polycarbonate 
cantilever in one of the two optical paths. The 
interferometer noise is about 5 x 10~ 14 m/VHz, 
which is two orders of magnitude smaller than the 
cantilever thermal fluctuations. The cantilever is 
inside an oven under vacuum. A window allows 
the laser beam to go inside (cf Fig.[2|). The size 
of the Polycarbonate cantilever are, I ~ 13.5 mm, 
a ~ 1 mm and b = 125 /jm, and the mirror mass is 

TOmirror ^$ 10 -3 g SUch that Pq ~ 100 Hz. 

Much care has been taken in order to isolate as 
much as possible the apparatus from the external 
mechanical and acoustic noise. The Nomarski in- 
terferometer and the cantilever are mounted on a 
plate which is suspended to a pendulum whose de- 
sign has been inspired by one of the isolating stages 
of the VIRGO superattenuator [HI E2 EH The 
whole ensemble is enclosed in a cage, to avoid any 
acoustic coupling. The pendulum and the cage are 




4Q Detector 



Oven under Vacuum Polycarbonate Cantilever (CI) 
plus its Glued Mirror 



Figure 2: Experimental setup. 



installed on air-suspended breadbord (Melles Griot 
Small Table Support System 07 OFA Series Ac- 
tive Isolation), which furnishes an extra isolating 
stage. However, these two isolation stages are not 
yet enough to have a large band measurement (0.1- 
100 Hz) of the cantilever thermal fluctuations. 

In Fig.Q] we plot the square root of the can- 
tilever's fluctuation spectral density as a function 
of frequency for a typical experiment at ambient 
temperature. The measure is compared with the 
FDT prediction, obtained from Eq.EH (blue line), 
and the interferometer noise (red line). The mea- 
sure scales quite well with the prediction. One can 
observe the cantilever resonance and the 1/y/v be- 
haviour for v <C vq produced by the viscoelastic 
damping (see Eci. I12J1 , However, the measurement 
is still too noisy in order to study accurately viola- 
tions of the FDT during aging. 

To improve our signal-to-noise ratio we have ap- 
plied our denoising technique described in Sec.|21 
As already mentioned, the total cantilever displace- 
ment reads x c = xt + x acc . To get xt we have to 
estimate x acc . The residual acceleration of the ta- 
ble where the interferometer is installed is about, 
10~ 8 m s~ 2 at 1 Hz and lO" 7 m s~ 2 at 100 Hz. This 
is too small to be detected by standard accelerom- 
eters, so we used a different method. We built an- 
other cantilever made by harmonic steel (cantilever 
C2) which is installed very close to the Polycar- 
bonate cantilever (cantilever CI). The parameters 




Figure 3: Square root of the fluctuation spectral pigure 4 . gquare root of the gpectral dengity of 
density without denoising (nm/^Hz) vs frequency Xaq yg frequency (Hz) . This signal is related to the 

*■ z '~ environmental noise. 



{a, b, I, m\ of C2 are chosen to optimize the sensitiv- 
ity to mechanical vibrations and reduce its sensitiv- 
ity to thermal noise (see Eo. lT5)l . A heavy mass and 
a heavy mirror, that give the main contribution to 
m for C2, are fixed on the steel cantilever free end. 
The cantilever C2 is damped by the viscosity of the 
air. A laser beam is reflected by the mirror glued 
on C2 and sent to a four quadrant position sensitive 
photodiode (4Q), which is used to detect the vibra- 
tions of the steel cantilever. The sensitivity is much 
smaller than that of the Nomarski interferometer, 
but enough for reducing the noise. Specifically, C2 
is 20 mm long, 10 mm wide, 0.125 mm thick, it 
has a total mass of 1.3 g approximately and a reso- 
nance frequency around 20 Hz. The maximum sen- 
sitivity external acceleration of this setup, which 
is limited by the 4Q detector, is about 10~ 7 m s -2 
in the frequency range of our interest. The out- 
put signal 1E4Q of the four quadrant detector and 
its Fourier transform X4Q are mainly proportional 
to the response of C2 to the external mechanical 
noise. Indeed, as we have already mentioned, ther- 
mal fluctuations of C2 are negligible. An example 
of the square root of the spectral density of the 4Q 
signal XiQ is plotted in Fig.^l which is related, via 
the response of C2, to the spectrum of the residual 
acceleration of the optical table. The polycarbon- 
ate cantilever and the steel cantilever, which are 



mounted very close on the same optical table, are 
perturbed by the same environmental noise sources. 
As theses sources may change of nature and of po- 
sition, the responses of CI and C2 to these external 
perturbations may change too. That is the reason 
why the denoising technique proposed in Sec.[2]can 
be very useful, because no hypothesis is needed on 
the response of the devices to the external noise. 
Referring to Sec.Gl one has to make the following 



substitutions: x 
whence 



X c , ^truo — * XT and £dc 



X±Q, 



(\x T \ 2 ) = (|i c | 2 )-|i? d | 2 (|i 4Q | 2 ), (16) 



with Rd 



L 4Q/ 



F4Q| 



(17) 



where the average (•) is computed in our experi- 
ment on a time interval T = 1 min, because Rd 
evolves on a time scale of a few minutes. This is 
shown in Fig.^lwhere we plot Rd, measured in three 
different time intervals separated by a few minutes. 
We see that the large variability of this response 
will make any a priori hypothesis useless. Using 
these data, we apply the denoising technique and 
we compute (\xt\ ) for each time interval of length 
T. Finally, we average the spectra obtained over 

several time intervals. 

In Fig.UOwe plot y (\xt\ ) as obtained after hav- 




Figure 5: Example of three |i? d | measured at dif- Fi S ure 6: Sc l uare root of the Ration spec- 

ferent time intervals of length T = 1 min separated tral denslt y wlth denoising (nm/ VHz) vs frequency 

by 2 min. Notice the large variation between the "> z '" 
three curves taken at different times. 



ing applied the noise reduction technique on twenty 
time intervals of length T = 1 min. By comparing 
this curve with Fig.Q]we see that all the peaks have 
been strongly reduced and that the agreement with 
the FDT prediction is much better than in Fig.^l 
Notice that no improvement is observed if the de- 
noising technique is applied on a single time inter- 
val of T = 20 min. We stress again that this effect 
is due to the fact that response is changing as a 
function of time. 

This example clearly shows that the denoising 
technique proposed in Sec.[21 can reduce the influ- 
ence of environmental noise on a measure if i?d is 
computed on short time intervals. The strong noise 
reduction introduced by this technique allows us to 
study the evolution of the FDT in an aging mate- 
rial. The accuracy is now limited by the 4Q noise, 
but this can be strongly reduced by replacing it 
with another Nomarski interferometer. 



5 Discussion and conclusions 

In this article, we have proposed an original and 
simple denoising technique, which allows one to re- 
duce the influence of the environmental noise on 
a measure. As already mentioned, this denoising 



technique mixes several features of the standard 
feedback systems and those of the Wiener filter- 
ing. The example presented in SecEl clearly shows 
that this technique can be very effective in sup- 
pressing spurious peaks on the spectra. The exam- 
ple of Sec.03 is not exhaustive. Indeed, the same 
technique can be used to reduce pick-up effects in 
electrical measurements or eventually in very pre- 
cise AFM measurements. 

As a general conclusions we can say that this 
technique is simple and can be implemented rather 
easily. The only requirement is to have a reli- 
able measurement of the environmental noise. Of 
course, it can be strongly improved by a multidi- 
rectional measurement of the noise. 
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